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1. INTRODUCTION 

Alternative current (AC) generators are considered as the main electrical energy sources particularly 
for high power systems. Power qualities such as stability, harmonics, and distortion are the main criteria of 
such energy sources. These qualities must be discussed particularly in induction generators which are 
recently considered as the main device in wind energy conversion systems due to their robust construction, 
low-cost maintenance and no separate excitation system [1]-[6]. From different factors having destructive 
effects on power qualities, magnetic saturation is considered as an important factor causing harmonic 
components through stator voltages and currents. It has been shown that the magnetic saturation produces 
destructive effects on performances of AC machines. The most effects are noise, increased losses and 
reduced insulation life [7]-[12]. In addition, the magnetic saturation causes the variation of machine 
inductances in terms of magnetizing current. 

Analysis of magnetic saturation in AC machines becomes an important topic of research and 
particularly attention is given to the modeling of magnetic saturation in induction machines. For this purpose, 
Roshandel et al. [13] investigated the use of subdomain technique to predict performance of the induction 
machine considering saturation effects. M6lsa et al. [14] proposed a dynamic model for saturated induction 
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machines with closed rotor slots and deep bars. Fatima ef al. [15] investigate the use of an improved 
permeance-based equivalent circuit model considering the frequency and magnetic curve dependent core loss 
branches for the prediction of the electromagnetic performances of induction machines. 
Mollaeian et al. [16] use fourier-based modeling of an induction machine considering the finite permeability 
and nonlinear magnetic properties. Ojaghi and Nashari [17] investigate the modeling of eccentric squirrel- 
cage induction motors with slotting effect and saturable teeth reluctances. Amiri et al. [18], Moreira and 
Lipo [19] investigate the modeling of saturated alternative current machines including air gap flux harmonic 
components. From several methods, the widely used is the finite elements method (FEM) [20]-[22]. This 
method is particularly characterized by the fact it gives good results, but its major inconvenience is the high 
computational time. The contribution of this work is that a state model has been developed for the healthy 
conditions of the squirrel cage induction generators (SCIG) considering the variation of machine inductances 
under the magnetic saturation. For this purpose, the machine inductances are calculated in terms of the 
magnetizing current. This approach is based on the assumption that both stator and rotor windings are 
replaced by an equivalent sinusoidally distributed windings [23], [24]. All the rotor variables are referred to 
the stator windings and called the equivalent variables. The rotor equivalent parameters are predicted in terms 
of the actual rotor bar resistances and inductances. 

To characterize the magnetic saturation effects in the SCIG, we are interested to the classical 
spectrum analysis due to its inherent advantages. This analysis is applied on stator currents, fluxes, and 
voltages. Experimental and theoretical results illustrate the consistency of this approach for the modeling and 
analysis of the magnetic saturation in the SCIG. 


2. SCIG EQUATIONS 
2.1. Stator voltage equations 

The direction of stator current is indicated by Figure 1. So, the stator voltage equation is expressed 
by (1). With the convention of the positive direction of stator and rotor currents, the stator flux becomes as 
given by (2). 


a[P,] _ ’ 

Pat = [Rll] + [m4] 
[®,] = [Pas Pps P5]" (1) 
[®,] = —[Lss] [is] = [Lsr] ili;] (2) 


The matrix of the stator inductances is expressed by (3). 


Lns this — bms/, > bms/, 
[Lss] aan iam bms/., Ims this = — bms/. (3) 
a bms/, = bms/, Lins + Lis 


Lms is the stator magnetizing inductance; and L,, is the stator leakage inductance. The magnetizing 
inductance is given by (4). 


2 
Lge be (4) 


Lm is the magnetizing inductance. Based on Kirchhoff laws, the stator and capacitor currents are expressed 
by (5) and (6) respectively. 


lisal = [ic] + [il (5) 
[i] = cS (6) 


For resistive load, the load current is given by (7). 
[i:] = [R)* [vn] (7) 
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C is the excitation capacitor. Z,, Z, and Z, are the load impedances. 
The matrix of load resistances is expressed by (8). 


rl O O 
[RI]J=|}0 rl O (8) 
0 0 rl 
So, the derivative of the vector of stator voltage is given by (9). 
d[vp : - 
Sal = = (Lig) — [Ril Iva) (9) 


[vn] = [Van Von Veal” 
lis] = [fas ips tal” 


Van> Vpn, and Ven are the phase voltages referred to the neutral point ‘n’. igs, i,,, and i,, are the line currents. 


lic] = twa ibe feel: 
Cis] = [igs ins tea)” 


lac, lye, and ic, are the currents in the capacitors. ig), ip;, and i; are the load currents. R; is the load 
resistance. 


Figure 1. Squirrel cage induction generator (SCIG) [24] 


2.2. Rotor voltage equations 

The squirrel cage rotor can be approximated as three identical windings. When expressing the rotor 
voltage equations, all rotor variables are referred to the stator windings and called the equivalent variables. 
With the convention of the positive direction of stator and rotor currents, the rotor voltage equation is 
expressed by (10). 


a[o,] _ 
— = [Rrlli-] (10) 


The rotor flux vector is expressed by (11). 


[,] = —[L;s] [is] = [Lye] [i] (11) 
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The matrix of the equivalent rotor-stator mutual inductances is defined by (12). 
[Ls] = [Lsr]” (12) 


[i,] is the vector of equivalent rotor currents; [R,-] is the matrix of the equivalent rotor resistances; [L,.,] is 
the matrix of the equivalent rotor-stator mutual inductances; [L,,] is the matrix of the equivalent rotor 
inductances. The matrices of the equivalent stator-rotor mutual inductances and rotor resistances are defined 
by (13) and (14) respectively. 


cos(@) cos(@ — =) cos(@ + =) 


[Ls] = Lms|cos(@ + =) cos(@) cos(@ — =) 


(13) 
cos(@ — =) cos(@ — =) cos(0) 
10 0 
[R,]=7,]0 1 0 (14) 
001 
For the matrix of equivalent rotor inductances, it’s expressed by (15). 
tims + Liy + bms/., _< bms/., 
= L L 
[Lr] a ia 5 Lins + Lirr = — ues (15) 
L L 
= nel 7 ms [5 Lins + Lir 


It can be demonstrated that the equivalent rotor resistance and leakage inductance can be calculated 
using (16)-(18) [23]. 


3 Nz 
y= ( =) eT [rb(1 — cos(a,)) + re] (16) 
_ 6 (\* 79 le 3 & 
Lin = 5 (5) N2 [2lo + srg] + 5Lms [Sag — 1] (i) 
6== (18) 


2.3. Calculation of machine currents 
To calculate the machine currents, it is convenient to determine the two components ias and ibs. The 
vector of the stator currents can be calculated using (19) to (21). 


[i,] = [Bs ][iaps] (19) 
With: 
+1 0O 
[B]=]0 41 (20) 
-1 -1 
tial = |," | 1) 


Using two independent components, we define the flux vector by (22). 


[Paps] = [As] [®,] (22) 


With the convention of the positive direction of stator and rotor currents, we can verify that the stator and 
rotor fluxes are calculated using (23) to (28). 
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[Pans] = —[Lsal liens] — Usral ir] (23) 

[%,] = —[Lysa] lias] — (Lr IL, ] (24) 
With, 

ieee Sl 0s) 

[Lea] = [As] [Ls] [Bsa] (26) 

[Lsra] = [As] [Lsr] (27) 

[L+sa] = [L,s][Bs] (28) 


Using (23) and (24), we get (29) and (30). 
Liaps] = —(Lsa] _ [Lsra] [Ler]* brsal)~* [Pans] = [Leva] [L,r]-1®,]) (29) 
[i,-] = —([L,,] “ [Lrsal [Lsal~*[Lsral)-*U®,] a [L+sallLsa]~*[Pans]) (30) 


2.4. SCIG state model 
Using (1), (9), (10), (29), and (30), the state model of SCIG is illustrated by (31). 


“es = —[Rs] [B,]([Lsa] eS [Lsra] Ler} Lrsa 7 * Pans] — [Lsra] [Lr] [8,]) ate [vy] 
AP] = TR, ]((Ler] — [Lysal bsal* seal) [Pr] — [rsallLsal*[®ans]) 
ee (31) 


dtc (IT] [Bs] [C5] Pans] ~~ Lsra] [Lr] 1) _ [Ril ttn) 
J = Tn — Te — f,a 


Where J is the inertia of the rotor and the connected load, 7. is the electromagnetic torque, 7, is the 
motorized torque, 2 the mechanical angular speed and fvis the viscose friction coefficient. The 
electromagnetic torque is expressed by (32) [25]. 


Pr. qt AlLsrl 


Te = 2 haiig 6 [i,] (32) 


P is the number of poles pairs and @ is the electrical angular displacement of the rotor. 


3. SIMULATION RESULTS 

Using the MATLAB environment, the above state model is simulated for a SCIG of 4 kW, 220/ 
380 V — 50 Hz. The excitation capacitor is C = 50 uF. A starting up with no load is considered. At time t = 
2.5 s, a balanced load of resistance 125 Q is supplied by the SCIG. The SCIG is rotated at mechanical speed 
of 3050 tr/mn imposing fundamental frequency of 50 Hz for the stator voltages and currents. The stator 
voltage and its spectrum analysis are illustrated by Figures 2(a) and 2(b), respectively. The stator flux, their 
hodograph and spectrum analysis are shown by Figures 2(c)-2(e), respectively. The spectrum analysis of the 
stator current is illustrated by Figure 2(f). 

As results of magnetic saturation effects, specific harmonics of order 2k+1 with k an integer 
(150 Hz, 250 Hz, 350 Hz, and 450 Hz...) are observed through the spectrum analysis of the stator voltages, 
fluxes, and currents as mentioned by Figures 2(b), 2(e), and 2(f). These harmonics explain the destructive 
effect of magnetic saturation such as noise, increase losses, reduced insulation life. As illustrated by 
Figure 2(a), Figure 2(c), and Figure 2(d), the temporary characteristics of the stator voltage, flux and the 
circular form of the stator flux are explained by the fact that the SCIG operate in balanced conditions. 
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Figure 2. Simulation results of SCIG: (a) stator voltage Van [24], (b) spectrum analysis of stator voltage Van, 
(c) stator flux phas [24], (d) hodograph of stator fluxes, (e) spectrum analysis of stator flux phas, 
and (f) spectrum analysis of stator current ias 


4. EXPERIMENTAL STUDY 
4.1. Magnetizing inductance 

To validate the theoretical study, an experimental test has been conducted on the three phase SCIG 
of 4 kW, 220/380 V — 50 Hz. The magnetizing inductance has been identified experimentally and expressed 
in terms of magnetizing current using an approximate function. For this purpose, the SCIG is coupled with a 
synchronous motor and a variable voltage is applied on the stator windings to measure the effective value of 
magnetizing flux in terms of the magnetizing current. Then, the magnetizing inductance is illustrated by 
Figure 3. In the literature, some functions have been proposed to approximate the magnetizing inductance in 
terms of the magnetizing current. In this study, it is approximated as indicated by (33). 

arctan(0.9im) 

Lm = ae (33) 
4.2. Validation of the theoretical study 

To validate the theoretical study, an experimental test has been conducted on the three phase SCIG 
of 4 kW, 220/380 V — 50 Hz. The excitation capacitor is C = 50 uF’. The SCIG is rotated at mechanical speed 
of 3050 tr/mn imposing fundamental frequency of 50 Hz for the stator voltages. A starting up with no load is 
considered and followed by the applying of balanced load of resistance 125 Q. The stator currents are 
measured using shunt resistors. The stator voltage and its spectrum analysis are illustrated by Figure 4(a) and 
Figure 4(b) respectively. 

It is important to note that the magnetizing inductance is sensitive to the magnetizing current 
depending on the load conditions. Figure 4(a) shows the temporary characteristic of the stator voltage which 
depends on the mechanical speed and the load conditions. Concerning Figure 4(b), it illustrates the presence 
of the harmonics (150 Hz, 250 Hz, 350 Hz...) through stator voltage. Systematically, these harmonics will be 
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reproduced through the stator currents and fluxes. These harmonic components will generate destructive 
effects such as noise, excessive heating and reduce insulation life of the SCIG. 
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Figure 4. Experimental results of SCIG: (a) stator voltage Van [24] and 
(b) spectrum analysis of stator voltage Van 
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5. CONCLUSION 

An analytical approach for the modeling and analysis of SCIG considering magnetic saturation has 
been synthetized and experimentally validated. The comparison of experimental and simulation results shows 
that harmonic components of the stator currents are produced due to the magnetic saturation. Using this 
approach, various characteristics are analyzed such as stator voltages and currents. This work must be 
performed by a future work to clearly discriminate the magnetic saturation from other stator asymmetries 
where all cases exhibit similar kind of current and flux signatures. 
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